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ABSTRACT 


A  summary  is  presented  of  progress  in  a  study  of  nonequilibrium 
generators,  with  emphasis  on  the  vibrational  excitation  of  molecular 
species,  experimental  vTork  on  a  large  generator  includes  development 
of  a  coaxial  preionizer  and  development  of  slanted  electrode  walls, 
both  of  which  increase  the  power  density  of  the  generator.  A  three- 
dimensional  analysis  of  the  slanted  electrode  wall  is  also  given. 
Finally  a  preliminary  report  is  given  of  experimental  and  theoretica1 
work  on  the  effect  of  an  intense  eradiation  field,  coupled  ho  the 
molecules,  on  electrothermal  instabilities. 
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I.  INTRODUCTION 


This  report  summarizes  progress  in  a  continuing  study  of  advanced 
concepts  in  energy  conversion.  The  general  objectives  of  the  program  are 
to  identify  promising  new  methods  of  energy  conversion,  and  to  carry  out 
the  research  required  to  determine  their  feasibility. 

The  work  has,  for  this  reporting  period,  concentrated  on  studies  of 
nonequilibrium  MHD  generators,  and  in  particular  on  two  aspects,  the  behavior 
of  slanted  electrode  walls,  and  the  behavior  of  molecular  species  in  the 
generator. 

As  noted  in  Annual  Report  No.  1  (l),the  slanted  wall  promises  to 
greatly  alleviate  the  problems  of  electrode  wall  shorting,  and  the  consequent 
limits  on  Hall  field  recovery.  Some  encouraging  preliminary  experimental 
results  were  given  in  Quarterly  Technical  Progress  Report  No.  22  (2).  A 
preliminary  finite-element  analysis  of  the  slant  wall  was  also  given  at  that 
time. 

Shortly  after  this  data  was  obtained  the  Inert  Gas  Heater  suffered  a 
mechanical  failure,  and  has  been  inoperable.  Repairs  were  begun  in  October 
1970,  and  it  will  be  in  operation  in  January  1971.  Meanwhile  analysis  of 
the  data,  and  more  detailed  theoretical  studies  of  the  slant  wall  have 
proceeded,  and  will  be  reported  here. 

In  a  parallel  pror .'am,  the  feasibility  of  vibrationally  exciting 
molecular  species  in  nonequilibrium  MHD  generators  wa3  demonstrated  (3). 

Because  of  the  Air  Force's  interest  in  such  possibilities,  the  generator 
program  was  reoriented  in  September  1970  to  place  emphasis  on  the  excitation 
of  molecular  species  in  nonequilibrium  generators.  Three  problems  were 
identified  for  study,  namely  1)  the  design  of  electrode  walls  suitable  for 
operation  at  the  desired  high  Hall  parameter  and  Mach  number,  2)  inlet 
relaxation  and  preionization  under  the  same  conditions,  3)  the  effect  on 
electrothermal  instabilities  of  strong  coupling,  through  molecules,  to  an 
intense  radiation  flux. 

The  first  two  items  are  essentially  continuations  of  past  work,  except 
that  emphasis  is  put  on  maximizing  incernai  dissipation,  rather  than  on 
producing  power.  The  last  is  n^w,  and  has  motivated  the  design  and  construction 
of  a  small  scale  experiment,  which  will  be  described.  Item  three  is  also 
being  studied  theoretically,  by  numerical  methods. 
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II.  SUMMARY  OF  PROGRESS 


2.1  Nonequilibrium  Generator  Experiment 

As  was  noted,  the  Inert  Gas  Heater  has  been  inoperable  since  May  1970. 

Our  effort  during  this  time  has  been  given  to  study  of  the  results  obtained 
prior  to  that  time,  to  conducting  a  cold  flow  test  on  the  channel,  and  to 
development  of  a  computer  program  for  calculation  of  the  flow  through  the 
channel.  The  following  results  have  emerged. 

Preionizer  behavior: 

The  preionizer  comprises  a  cathode  near  the  throat  of  the  axially  seg¬ 
mented  supersonic  nozzle,  and  an  anode  near  the  channel  entrance.  A  "spinning 
coil"  is  arranged  to  impress  an  axial  magnetic  field  on  the  nozzle.  It  vas 
observed  that,  with  this  spinning  coil  in  operation,  the  channel  produced 
near-ideal  open  circuit  voltage,  but  when  loaded  it  performed  very  poorly. 
Movies  taken  from  the  channel  exit  showed  that  this  was  due  to  a  stratification 
of  the  plasma,  in  the  UXB  direction,  produced  by  the  preionizer  current 
flowing  along  magnetic  field  lines  from  the  cathode  tip  to  the  insulator  walls. 
A  frame  of  this  movie,  shown  in  Fig.  (la)  exhibits  this  pattern.  Fig.  (lb) 
shows  the  same  view  with  the  spinning  coil  off.  It  exhibits  a  more  uniform 
plasma,  because  in  this  operating  mode  the  discharge  was  from  the  cathode  tip 
to  the  nozzle  throat.  We  conclude  that  the  stratification  reduced  the  effect 
conductivity  in  the  channel,  preventing  large  load  currents.  The  preionizer 
has  been  modified  by  addition  of  anode  ballast  resistors  to  ensure  a  more 
uniform  current  distribution. 

Gas  dynamics: 

There  was  some  uncertainty  in  the  gas  dynamic  behavior  of  the  channel, 
with  considerable  evidence  that  shocks  existed  in  the  downstream  section  (2). 
Accordingly  the  channel  side  walls  were  replaced  with  plexiglass  plates ,  and 
the  complete  assembly  of  preionizer,  channel,  and  diffuser  was  tested  in  a 
supersonic  wind  tunnel  loop.  The  channel  is  shown  in  place  in  Fig.  (2.).  A 
shadowgraph  of  the  upstream  portion  of  the  channel.  Fig.  (3a),  shows  a  weak 
oblique  shock  system  originating  from  the  preionizer  nozzle.  This  was  ex¬ 
pected.  Fig.  (3b)  shows  two  normal  shocks  in  the  downstream  portion  of  the 
channel.  Their  locations  agree  with  those  inferred  from  wall  pressure 
measurements.  This  offers  an  explanation  of  abrupt  drops  in  open-circuit 
voltage  reported  in  Ref.  (1),  as  well  as  a  cutoff  of  load  current  in  the 
downstream  half  of  the  channel. 

Electrode  drops: 

Comparison  of  the  data  from  the  flat  wall  and  slant  wall  channels  shows 
electrode  drops  of  approximately  10  volts  for  the  slant  wall  compared  to 
approximately  1*0  volts  for  the  flat  wall.  This  improvement  evidently  results 
from  current  flow  along  the  magnetic  field  lines. 

2.2  Three  Dimensional  Analysis  of  Slanted  Electrode  Wall 

The  basic  idea  of  the  slanted  elect'-ode  wall  is  to  relieve  tne  boundary 
condition  imposed  by  the  insulating  segments  of  the  electrode  wall,  by 
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allowing  a  flow  of  current  to  the  wall,  in  the  direction  of  the  magnetic 
field.  This  is  a  three  dimensional  phenomenon,  yet  analyses  of  it  have 
been  essentially  two-dimensional.  A  three  dimensional  analysis  of  the 
current  and  field  structure  is  given  in  detail  in  section  3.1,  for  the 
case  of  constant  conductivity.  The  analysis  shows  that  a  condition  for 
minimum  dissipation  on  the  electrode  wall  is  that  the  field  along  the  in¬ 
sulating  segments  in  the  direction  perpendicular  to  the  generator  axis 
should  be  minimized.  This  requires  narrow  insulators  and  (or)  fine  seg¬ 
mentation.  The  calculations  show  that  the  slant  wall  does  in  faet  reduce 
the  current  near  the  electrode  leading  edge,  as  expected.  Strong  (rotation) 
eddy  currents  are  predicted  near  the  apex  of  the  slant  wall,  which  might 
lead  to  substantial  dissipation  in  the  nonequilibrium  plasma.  Alternate 
wall  shapes  are  being  sought  to  minimize  these  effects. 

2.3  Experimental  Study  of  the  Interaction  of  Electrothermal  Instabilities 

with  an  Insulated  Obstacle 

An  experiment  was  conducted  to  simulate  the  behavior  near  a  segmented- 
electrode  wall,  by  inserting  an  insulating  obstacle  in  a  square  channel  with 
axial  current  flew  and  magnetic  field  normal  to  its  axis.  The  spaces  between 
the  obstacle  and  the  tube  walls  simulate  the  electrode  surfaces,  as  current 
passes  through  them,  while  the  obstacle  itself  simulates  the  Insulator  seg¬ 
ment.  A  computer  simulation  of  the  same  geometry  was  also  carried  out,  with 
the  aim  of  direct  comparison. 

The  results  of  the  experiment  are  confusing.  Apparently  the  (unsegmented) 
cathode  and  anode  at  the  ends  of  tube  influenced  the  current  pattern  more 
than  was  expected,  so  that  direct  comparison  of  the  experiment  and  calculation 
has  not  so  far  been  possible.  Both  aspects  will  be  described  in  detail  in 
the  next  quarterly  report. 

2.4  Pulsed  Experiment  with  Molecular  Gases 

In  order  to  determine  the  effect  of  an  intense  radiation  field  on  the 
behavior  of  a  nonequilibrium  plasma  with  molecular  species,  a  small,  short- 
time  scale  experiment  has  been  designed.  It  consists  of  a  quartz  discharge 
tube,  shown  in  Fig.  (4),  with  end  mirrors,  a  pulsed  air  core  magnet,  and 
apparatus  for  supplying  the  tube  with  He,  Cs,  CO  or  CO2  mixtures.  The  mode 
of  operation  will  be  as  follows:  A  slow  flow  of  seeded  He-C02  mixture  will 
be  passed  through  the  tube  as  indicated  in  the  flow  schematic  of  Fig.  (5). 

An  axial  condenser  discharge  will  produce  the  plasma  on  a  time  scale  long 
compared  to  that  required  to  produce  the  ionization  and  to  excite  vibrational 
states,  but  short  compared  to  that  required  to  heat  the  gas.  A  schematic  of 
the  electrical  system  is  shown  in  Fig.  (6).  Measurements  will  include  axial 
and  Hall  fields,  Cs  continuum  radiation,  and  axial  10. 6p  CO2  radiation.  Hie 
major  components  of  this  apparatus  are  shown  in  Fig.  (7).  Experiments  will 
begin  during  February  or  March,  1971. 

2.5  Ionization  Instabilities  in  a  Magnetohydrodynamic  Laser 

The  two  temperature  operating  regime  of  the  MHD  laser  will  be  one  in 
which  'onization  instabilities  would  normally  be  present.  These  ionization 
instabilities  will  be  affected  by  the  presence  of  the  lasing  radiation  field. 


In  a  two  temperature  plasma  with  molecular  species  but  without  a  radiation 
field  which  interacts  with  the  molecules,  ionization  instabilities  are  not 
significaily  modified  by  the  molecular  species.  If  the  excited  states  of 
the  molecular  species  are  in  equilibrium  at  the  gas  temperature,  the 
molecular  species  simply  increase  the  electron  heavy  particle  loss  rate 
without  qualitatively  changing  the  nature  of  the  instability.  If  the 
molecular  species  excited  states  are  in  equilibrium  at  the  electron  tem¬ 
perature,  then  the  effect  on  the  electron  gas  may  be  described  in  terms  of 
a  somewhat  larger  specific  heat.  The  onset  condition  of  the  instability 
should  not  change  for  this  case  but  the  growth  time  is  lengthened  slightly. 

When  the  molecular  species  interacts  with  a  lasing  radiation  field, 
one  expects  qualitative  changes  in  the  behavior  of  the  instability.  This 
is  because  the  electron  gas  is  now  coupled  through  the  molecular  species 
to  the  radiation  field  and  a  spatial  redistribution  of  the  electron  energy 
is  possible.  Further,  if  ionization  instabilities  persist  in  the  presence 
of  the  radiation  field,  the  radiation  field  will  be  subject  to  a  'turbulent' 
gain  dist^ioution  within  the  plasma  stemming  directly  from  the  nonuniformities 
resulting  from  the  ionization  instability. 

The  appropriate  rate  and  conservation  equation  for  the  CO2  vibrational 
system  and  the  electron  energy  equation  including  electron  excitation  of  the 
CO2  lasing  vibrational  levels  have  been  formulated  along  with  the  Maxwell 
equations  describing  the  electrical  conduction  field  in  the  MHD  laser  and 
the  radiative  transport  equation  describing  the  radiation  field.  This 
description  will  now  be  used  to  study  both  one  and  multi-dimensional 
aspects  of  the  coupling  of  the  two  temperature  plasma  and  the  lasing 
radiation  field.  The  effect  of  the  radiation  field  on  the  ionization  in¬ 
stability  will  then  be  examined  as  well  as  the  effect  of  the  turbulent 
plasma  on  the  radiation  field. 
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III.  DETAILED  TECHNICAL  CONTRIBUTIONS 


3.1  The  Three-Dimensional  .electrical  Conduction  Field  in  a  Slanted 

Electrode  Wall  MHD  Generator 

3.1.1.  Introduction 

In  a  conventional  finitely  segmented  Faraday  connecteu  MHD  generator, 
current  concentration  effects  at  electrode  edges  and  ov*r  insulators  exist 
vhich  lead  to  progressively  stronger  Joule  dissipation  . .  these  regions  as  the 
Hall  parameter  is  increased.  These  geometrical  nonuni fortuities ,  induced  by 
the  adjacent  locations  of  insulator  and  conductor  segments  were  first  pre¬ 
dicted  by  Hurvitz,  Sutton,  and  Kilb  (h),  and  have  been  a  continual  under¬ 
lying  limitation  to  the  electrode  wall  performance  of  MHD  channels.  Although 
progressively  finer  segmentation  may  alleviate  the  deleterious  effects  of 
the  conductor  insulator  nonuniformity  in  equilibrium  gas  generators  operating 
with  modest  Hall  parameters,  Kerrebrock  (5)  has  shown  that  these  nonuniformities 
exert  a  more  serious  effect  on  nonequilibrium  generators  where  the  excess 
Joule  dissipation  at  the  electrode  wall  readily  elevates  the  conductivity  of 
the  gas  near  the  wall  thereby  rendering  the  wall  segmentation  ineffective. 

Ricateau  (6)  has  suggested  that  the  current  concentrations  in  MHD 
channels  with  electrode  walls  parallel  to  the  magnetic  field  may  be  alleviated 
by  inclining  these  electrode  walls  to  the  magnetic  field.  In  this  way,  the 
current  is  free  to  enter  the  electrode  along  the  magnetic  field  lines  and  need 
not  concentrate  at  the  electrode  insulator  intersections  as  in  the  parallel 
wall  case.  In  the  present  work  we  present  a  detailed  study  of  the  electrical 
conduction  field  in  such  a  slanted  wall  channel.  In  Part  3.1.2  we  examine 
the  conditions  for  minimum  Joule  dissipation  on  insulating  and  conducting 
surfaces  inclined  to  the  magnetic  field.  In  Part  3.1.3  we  present  the  de¬ 
tailed  three-dimensional  structure  of  the  fields  and  currents  in  a  slanted 
wall  MHD  channel.  In  Part  3.1.^  we  examine  the  effect  of  thick  velocity 
boundary  layers  on  the  behavior  of  the  current  in  the  slanted  wall  region. 

In  Part  3.1.5  we  show  that  Hall -effect  induced  eddy  currents  resulting  from 
nonuniformities  in  the  magnetic  field  direction  exist  even  for  the  case  of  a 
uniform  gas. 


3.1.2  Joule  Dissipation  on  Surfaces  Inclined  to  the  Magnetic  Field 

We  consider  a^gas  with  conductivity*”?  supporting  a  current  ?  with 
induced  field  E'  =  E  +  IJ  x  $  where  for  Cartesian  coordinates  x,  y,  z  with  B 
parallel  to  the  z  axis , 

?  »  o?  E’  (1) 


The  conductivity  tensor  is  given  by 


o  = 


-So, 


1  +  6 

So* 

1  +  3* 

0 


1  +  6 

-A— 
1  +  S2 

0 


(2) 
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where  a  and  o  are  conductivity  components  and  8  is  the  Hall  parameter. 

If  ve  transform  from  current  components  (Jx,  Jy ,  Jz)  to  components  (Jx,  Jn, 
Jt)  where  n  and  t  denote  direction  normal  and  tangential  to  an  electrode 
wall  surface  inclined  to  the  8  field  at  angle  a,  ve  obtain 


-Bcosa 


-Bsina 


J  =  — Bcosa  cos2a+-^(l+B2)sin2a  [l — ^(l+82)]sinacosa  E' 

n  1+02  °x  °x  n 

J.  Bsina  [l---^(l+82)]£inacosa  sin2a+-^(l+B2)cosa  E' 


On  surfaces  where  Jn=0  (insulating  surfaces),  we  find  from  (6)  that  the 
normal  field  and  transverse  and  axial  currents  are 

2 

_  Bcosa  ,  B  sinacosa  , 

n  2  .  2  Ex  +  1a02  “2  Et 

1+B  sin  a  1+8  sin 


T  s  -  -  u 

t  _02  .  2 
1+8  sin  a 


( BsinaE'+E' ) 

X  T* 


2  2 
1+8  sin  a 


(E'  ~8sin  EM 

X  t 


where  we  have  assumed  0|,=Oi,=o.  It  is  noteworthy  that  the  last  two  of  Eqs. 
(4)  are  formally  identical  to  the  usual  Ohm's  Law  in  the  plane  perpendicular 
to  B  with  8  replaced  by  Bsina,  i.e.,  the  Hall  parameter  is  based  upon  the 
component  of  8  normal  to  the  plane  in  which  the  current  vector  lies.  This 
formal  similarity  provides  the  clue  to  a  requirement  for  minimum  Joule 
dissipation  on  insulating  surfaces  inclined  through  the  magnetic  field: 
given  a  field  E'  impressed  on  the  surface  by  the  core  of  the  channel,  it  is 
necessary  that  $he  "Hall  field"  E£  be  made  to  vanish.  This  requirement  is 
clearly  seen  when  the  Joule  dissipation  on  the  insulating  surfaces  is  con¬ 
sidered: 

2  2 

.2  E'  +E? 

J  x  t 

s  n  ■— 

<7  -2  .2  / 

1+8  sin  a  (5) 


Thus  contrary  to  the  usual  situation  in  the  core  of  the  channel  where 
maximum  recovery  of  Hall  field  is  desirable,  on  the  inclined  insulating 
wall,  it  is  desirable  that  the  appropriate  Hall  field  E£  not  be  recovered 
at  all.  The  elimination  of  the  Hall  field  E{  may  be  achieved  by  arranging 
conducting  and  insulating  segments  with  fine  segmentation  relative  to  the 
channel  width  so  that  the  close  proximity  of  conductors  imposes  a  tangential 
equipotential  on  adjacent  insulating  surfaces.  We  note  that  if  E{=0,  a 
tangential  current  will  flow  on  the  insulator  surface  given  by 


6 


wwi»iw 


J 


t 


gsina 

lx.2  .  2 

1+6  sin  a 


a  E’ 
x 


(6) 


If  the  imposition  of  the  tangential  equipotential  is  ineffective  and  the 
tangential  current  J+  is  impeded  so  that  J+s0,  the  tangential  field  which 
develops  is  E£  =  -gsinaE^.  The  Joule  dissipation  in  this  case  is 


-  l+6^sin^o(— ) 
j  *o  0  e;=o 

t  t 


2  2 

and  for  large  6  is  of  order  of  6  sin  a  times  the  dissipation  with  E^=0. 

If  the  insulating  wall  is  in  contact  with  a  core  in  which  full  Hall  recovery 
exists  and  Ex=(Ex)core=inJy)core,  t^en  t^e  dissipation  on  the  insulator 
wall  relative  to  the  coFe  Is 


2  2 

1+6  sin  a 


core 


(~)  =  e2  (~) 

a  T  n  o 

J. =0  core 

”C 


(7) 


Equations  (7)  exhibit  the  crucial  features  of  slanted  wall  elec¬ 
trodes.  For  6  »  1,  the  slanted  wall  (a^O)  leads  to  an  insulator  wall 
dissipation  which  is  of  the  order  of  the  core  dissipation  (provided  E£=0) 
whereas  in  the  parallel  wall  case  (a=0),  the  insulator  wall  dissipation  is 
of  the  order  6^  times  the  core  dissipation.  If  the  condition  E{-0  is  not 
secured  and  in  fact  the  current  is  impeded,  the  effect  of  the  slanted 
wall  in  diminishing  the  insulator  wall  dissipation  to  that  of  the  core  will 
be  lost  even  if  a^O. 


The  situation  on  a  conducting  surface  may  be  similarly  examined. 
Imposing  the  equipotential  conditions  Ex  =  0,  =  0  and  assuming  E-t  =  E£ 

near  the  wall  (the  thick  boundary  layer  case),  the  normal  and  axial  currents 
are  found  from  Eq.  (3)  as 


J 
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-6‘^sinacosa 

2  2 
1+6  sm  a 

-6  cosa 

*  2  2 

1+6  sm  a 


(8) 


The  Joule 


o 


dissipation  on  the  conducting  surface  is  then 


1+6 


2  .  2 
1+6  sin  a 


If  it  is  again  assumed  that  the  conducting  surface  is  in  contact  with  a  core 
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with  full  recovery  of  the  Hall  voltage  and  that  J  =(J  )  ,  we  find  that 

n  y  core 

2  2  ? 

J_  =  1+B  ( J^j 

0  l+0^sin^a  0  core  (9) 

The  conducting  surface  in  a  slanted  wall  channel  will  therefore  exhibit 
dissipation  on  the  order  of  the  core  dissipation  whereas  with  o  =  0,  the 
conductor  surface  dissipation  is  of  the  order  B2  times  the  core  dissipation 
for  B  »  1. 


3.1.3  Detailed  Structure  of  the  Electrical  Conduction  Field 


In  a  finite  sized  channel  the  behavior  of  the  current  on  the  con¬ 
ducting  and  insulating  surfaces  must  be  completed  by  consideration  of  the 
nonuniform  boundary  conditions  associated  with  the  location  of  adjacent 
insulators  and  conductors  and  by  the  presence  of  insulating  sidewalls  in  the 
magnetic  field  direction.  We  shall  therefore  consider  the  structure  of  the 
three-dimensional  electrical  conduction  field  subject  to  the  boundary  con¬ 
ditions  imposed  by  a  finite  sized  channel  with  finite  segmentation.  The 
Maxwell  equations  for  the  small  magnetic  Reynolds  number  regime  and  a  time 
scale  for  conduction  long  compared  to  the  plasma  frequency  period  are 


Vjli  a  0 


v  -  0 


(10) 


We  assume  £he  velocity  field  is  such  that  V  x  2'=  0  and  potentials  4>  exist 
such  that  E'  ■  -V*',  E  ■  -V*.  For  the  case  of  uniform  a,  8  Eqs.  (l),  (3) 
and  (10)  reduce  to  a  single  elliptic  equation  for  the  potential 


3-f  ♦  ♦.*)*% 

3x  3x  3z 


We  shall  apply  Eq.  (11 )  to  the  description  of  the  fields  and  currents  which 
exist  in  the  chevron  shaped  slanted  wall  channel  structure  of  Fig.  1.  The 
boundary  conditions  applied  are  those  of  vanishing  normal  current  on  insula¬ 
ting  surfaces  and  specified  potentials  on  conducting  surfaces.  A  single 
section  consisting  of  one  electrode  pair  of  an  infinitely  long  periodic 
channel  may  be  described  by  four  global  variables:  the  potential  difference 
between  the  upper  and  lower  electrodes  V  ,  the  potential  difference  developed 
in  one  period  in  the  axial  direction  Vx,^the  total  current  flowing  in  the 
electrode  pair  Iy  and  the  total  current  flowing  in  the  channel  in  the  axial 
direction  Ix.  As  shown  in  Ref.  (9),  any  suitable  pair  of  these  global 
variables  Vx,  Vy ,  Ix,  Ty  determine  the  complete  distribution  of  field  and 
current  within  the  channel. 


Ou  insulating  surfaces  perpendicular  to  the  magnetic  field,  the 
appropriate  boundary  condition  is  3$J_  .  On  the  insulating  surface  within 

dz 
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the  slanted  electrode  vail,  the  condition  for  vanishing  normal  current  is 


.3*'  3*'  .  „ 

B3T  +  37  +  tana  (l+B  )a7  =  0 


The  periodicity  condition  is 


*'(x,y,z)  =  «'(x+t,y,z)-Vx 


(12) 


(13) 


The  upper  conducting  surface  is  at  a  constant  potential  <t>  = 
The  potentials  <t>  are  related  (for  uniform  magnetic  field)  by 


-V  /  2. 

y 


4>'(x,y,z)  -  4>(x,y,z)  +  B£yu(n)dn 


(HO 


where  0'(x,0,z)  =  4>(x,0,z) .  It  follows,  therefore,  that  on  the  conducting 
surface,  the  potential  4>'  is  given  by 


♦  ' 


=  -V  /2  +  B/yifii)dn 


(15) 


where  y(z)  is  the  y  coordinate  of  the  slanted  wall  given  by  y(z)=H/2+(a/2-z) 
tern  a.  The  potential  distribution  for  a  given  velocity  distrioution  may 
be  obtained  by  substituting  the  given  velocity  profile  into  Eq.  (15).  We 
shall  consider  only  two  general  cases:  that  of  thin  boundary  layers  relative 
to  the  channel  thickness  a  and  that  of  thick  boundary  layers  relative  to  a. 
For  the  former  we  set  u(n)  =  U  and  thus  determine 


*'  =  -V  /2  +  UBH/2  +  uB  (a/2  -  z)tancx 

y 


(16) 


For  the  latter,  we  set 


U 


0  ,<  y  {  II 


u(n)  = 


,0  H<  y 

so  that  in  the  thick  boundary  layer  case 


=  -V  /2  +  UBH/2 

y 

In  terms  of  the  potential  4>'  the  global  currents  I  ,T  are 

x  y 


(1Y) 
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I*  "  ~S  ' 
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o  ,3#'  3*\,  ,, 

T  (3T  '  8  sT)dsr  dz 


a/2  1/2 


y  -a/2  1/2  1+82  '  3x  8y 


o  3*'  3*' »  . 
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and  the  Joule  dissipation  is 


J2  _  „r  1  r/3*\2.  ,3*\2,  .  ,3*\2, 

a{  [(_)  ♦(—)]♦  (~)  } 


(18) 


(19) 


Since  the  problem  posed  is  linear  for  a  homogeneous  differential  equation 
with  partially  inhomogeneous  boundary  conditions,  the  general  solution  for 
any  pair  of  Vx,Vy,Ix,Iy  may  be  obtained  by  superposition  of  a  fundamental 
set.  This  fundamental  set  of  solutions  is  obtained  numerically  by  iteratively 
solving  a  set  of  finite  difference  approximations  to  Eq.  (ll)  and  a  fundamen¬ 
tal  set  of  boundary  conditions  based  upon  Eqs.  (12),  (13),  and  (l6).  Details 
of  the  numerical  technique  may  be  found  in  Ref.  (7). 

We  shall  describe  here  only  results  for  a  Faraday  connected  channel 
(Ix=0).  For  such  a  channel  two  performance  parameters  are  of  particular 
interest:  the  internal  impedance  R^  and  the  apparent  Hall  parameter  8app. 
These  are  defined  as 


app 


<E  >/<E’>  =  (V  /V  )(H/0 
x  y  x  y 


where  <  >  denotes  that  each  component  of  the  electric  field  is  averaged  in 
its  coordinate  direction.  For  a  three-dimensional  channel,  the  fields  must 
also  be  averaged  in  the  magnetic  field  direction.  For  a  simple  arithmetic 
average,  the  two  expressions  for  Sapp  corresponding  to  thick  and  thin  boundary 
layers  are: 


thin: 


app 


V 

x 


uBH-V  --  ln(l+a/H  tana)/tana 


Vx(a/il)  tana 

thick:  eapp  =  (-u3h_v  ) ( i/V)"3^1+a/H“ tanoT 

y 


The  results  to  be  described  are  for  a  modest  Hall  parameter  3=U  and  for  a 
slant  wall  angle  a=!45°. 
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Let  us  consider  first  the  behavior  of  a  slanted  va^  c  uel  at 
open  circuit.  In  the  thin  boundary  layer  case  the  induced  field  E'  will 
be  shorted  by  the  presence  of  the  conductors.  The  shorting  current  appears 
as  a  large  loop  current  which  leaves  the  anode  near  the  apex  region  and 
enters  near  the  sidewall.  For  the  cathode  the  pattern  is  reversed.  To 
study  the  loop  current,  the  open  circuit  problem  was  solved  for  a  fine 
electrode  segmentation  (£/H=l/10).  The  structure  of  the  open  circuit  loop 
current  is  shown  in  Fig.  6.  The  y-z  cross  section  shown  intersects  the 
electrode.  The  current  vectors  associated  with  the  loop  current  which  lie 
in  the  plane  of  the  electrode  wall  are  shown  in  Fig.  7.  The  large  positive 
x  components  of  tne  vectors  indicates  that  the  loop  is  highly  three-dimensional. 
Since  the  net  axial  current  Ix  must  be  zero,  the  loop  reverses  its  x 
direction  at  some  depth  within  the  channel .  The  distribution  of  normal 
current  on  the  electrode  wall  is  shown  in  Fig.  8.  Since  this  internal 
circulating  loop  current  remains  the  same  for  all  loads,  it  becomes  a  small 
fraction  of  the  total  generator  current  at  moderate  loads  if  the  channel 
thickness  a  is  small  relative  to  the  channel  height  H. 


It  is  convenient  to  describe  the  loaded  channel  in  terms  of  a 
load  factor  ky  =  Vy/UBH.  When  the  channel  is  under  load,  the  load  current 
dominates  the  loop  current  and  the  normal  current  distribution  shown  in 
Fig.  9  is  obtained.  It  can  be  seen  that  although  current  concentration? 
at  the  electrode  edges  exist,  these  concentrations  are  considerably  weaker 
than  tnose  in  parallel  wall  channels  at  the  equivalent  segmentation  and  Hall 
parameter.  It  may  also  be  noted  that  on  the  edge  with  the  weak  singularity, 
current  continues  to  leave  the  conductor.  This  is  a  result  of  both  the 
circulating  current  and  the  presence  of  "leakage"  [Crown  (8),  Oliver  (9)]. 
The  Joule  dissipation  shown  in  Fig.  9  is  similarly  dimished  over  that  for 
parallel  wall  channels  except  at  the  apex  where  strong  dissipation  exists. 
This  strong  dissipation  in  the  apex  region  is  due  to  Hall  effect  induced 
eddy  currents  which  are  circulating  in  the  apex  region.  These  eddy  currents 
will  be  discussed  in  Part  3.1.5.  It  cam  be  seen  that  the  potential  <J>  is 
quite  uniform  over  the  insulator  region  in  the  tangential  direction  running 
from  the  apex  to  the  side  wall.  The  field  is  therefore  shorted  out  and, 
as  discussed  in  Part  j.1.2,  minimum  dissipation  on  the  insulator  wall  is 
closely  achieved.  This  uniform  potential  condition  results  from  the  rather 
fine  segmentation  relative  to  the  channel  width  which  allows  the  adjacent 
conductor  to  impose  a  tangential  equipotential  on  the  insulator.  These 
minimum  dissipation  effects  are  reflected  in  the  performance  characteristics 


of  this  channel: 


3  =  3.2U,  R./(R. ). ,  .  =  1.79 

app  ii  ideal 


where  (R.  ),  ,  ,  is  the  internal  impedance  of  an  infinitely  fine  segmented 

,  ,  l  ideal 

channel . 


The  current  vectors  on  the  slanted  wall  are  shown  in  Fig.  10. 

It  should  be  noted  that  the  x  component  of  these  vectors  is  small  indicating 
the  effectiveness  of  the  slanted  wall  in  inhibiting  the  axial  Hall  current. 
The  current  patterns  in  the  channel  are  displayed  for  several  y-z  cross 
sections  in  Fig.  11.  Two  trends  are  present  depending  on  whether  the  cross 
section  passes  through  the  electrode  or  through  the  insulator.  For  the 
electrode  cross  section  the  current  leaves  the  electrode  with  a  negative  z 


component.  As  the  current  traverses  the  channel,  the  current  lines  curve 
until  the  z  component  is  zero.  Fcr  the  insulator  cross  sections  the 
current  initially  flows  along  the  insulator  vith  a  positive  z  component. 

The  eunent  lines  again  curve  until  the  z  component  is  zero.  In  both 
instances  the  curvature  occurs  within  a  distance  of  the  order  of  the 
electrode  width. 

A  three-dimensional  presentation  of  the  current  vectors  within 
the  channel  is  shown  in  Fig.  1?.  The  cross  sections  are  in  the  x-y  plane 
at  four  z  locations.  The  sections  extend  only  a  depth  l  into  the  channel. 
The  variation  in  the  z  component  of  the  current  vectors  evident  in  Fig.  11 
is  also  present  here.  An  examination  of  the  rows  of  vectors  beneath  the 
"hot  edge"  shows  how  sharp  the  reversal  of  the  z  component  is. 

When  these  same  calculations  are  carried  out  for  a  more  coarsely 
segmented  channel  (t/a  -  1/2),  it  is  found  that  the  internal  impedance 
increases  to  Ri/(Ri) ideal  =  5.22  while  the  apparent  Hall  parameter  falls 
to  3apps2.0U.  The  potential  and  current  on  the  electrode  wall  are  shown 
in  Fig.  13  for  this  case.  It  can  be  seen  that  the  Joule  heating  is  more 
intense  and  the  tangential  equipotential  on  the  insulator  wall  fails 
consistent  with  the  deterioration  in  channel  performance.  The  current 
vector  patterns  on  the  slanted  wall  (Fig.  lU)  are  similar  to  those  for 
the  fine  segmentation  case;  however  the  vectors  for  the  coarse  segmentation 
case  have  larger  x  components  associated  with  them.  The  current  patterns 
in  the  y-z  cross  sections  re  shown  in  Fig.  15.  The  distance  over  which 
the  current  vectors  turn  to  become  perpendicular  to  the  magnetic  field  (z) 
direction  is  larger  than  in  the  fine  segmentation  case  but  still  of  the 
order  of  the  electroce  width. 

3.1.4  Effect  of  Thick  Boundary  Layers 

If  the  velocity  boundary  layers  are  thick  relative  to  the  channel 
thickness  so  that  the  fluid  is  essentially  stagnant  in  the  slanted  wall 
region,  then  for  such  a  situation,  the  open  circuit  loop  current  will  not 
exist  and  the  shorting  effect  of  E'  even  under  load  will  not  be  present. 
Calculations  were  carried  out  for  the  fine  segmentation  case  (t/a=l/10). 

The  values  of  internal  impedance  and  apparent  Hall  parameter  for  the  thick 
layer  case  were  3app=3.20  and  Pi/(Ri)ideal  =  1.19.  Although  the  apparent 
Hall  parameter  is  little  changed  from  the  thin  boundary  layer  case,  the 
elimination  of  the  internal  circulating  current  improves  the  internal  im¬ 
pedance  significantly  for  the  considered  geometry  in  which  a  =  H.  The 
normal  current  distribution  in  Fig.  16  exhibits  "leakage"  indicating  that 
the  simultaneous  emitter  and  collector  characteristics  of  the  electrode 
were  not  due  only  to  the  loop  current.  Analogous  to  the  two-dimensional 
stagnation  point  is  the  line  of  points  with  zero  current  shown  in  Fig.  7. 

3.1.5  The  Apex  Region 

The  apex  region  for  the  chevron  shaped  channel  is  a  singular 
region  where  a  sharp  boimdary  discontinuity  in  the  z  uirection  occurs, 
because  of  the  presence  of  the  central  symmetry  plane  z  =  0  across  whicn 
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there  is  no  current  flow,  the  normal  and  tangential  current  components 
diminish  to  very  small  values  near  the  apex;  however,  because  of  the  Hall 
effect,  strong  axial  current  may  still  flow  in  the  apex  region.  To  improve 
the  resolution  near  the  apex,  a  finer  grid  was  employed  in  that  region. 

The  results  of  these  refined  calculations  are  shown  in  Fig.  19.  It  can  be 
seen  that  the  peaking  and  diminishment  of  the  normal  current  as  the  apex 
is  approached  is  revealed.  The  Joule  dissipation  however  exhibits  strong 
singularities  on  the  insulator.  These  intense  regions  of  Joule  heating 
result  from  eddy  currents  which  are  circulating  in  the  apex  region. 


The  circulation  of  the  current  field  may  be  determined  from  Eqs. 
(l)  and  (10)  as 


V  J  =  -B*vJ  +  J*V0  +  (J-Vlno)B  -  [(3 ‘Vino  -  V*0]j 


(20) 


In  the  case  of  a  general  nonuniform  gas,  the  circulation  of  the  current 
field  is  proportional  to  gradients  in  the  conductivity  components;  however, 
even  for  a  uniform  gas,  the  presence  of  boundary  induced  nonuni  formitie*? 
in  the  direction  of  the  magnetic  field  will  lead  to  a  rotational  current 
field  given  by  the  first  term  on  the  right  hand  side  of  Eq.  (20).  We  note 
however  that  according  to  Eq.  (20)  a  uniform  gas  conducting  current  only 
in  the  plane  perpendicular  to  0  will  be  free  of  circulation.  In  Figs.  20 
and  21  the  current  vectors  in  the  plane  of  symmetry  z  *  0  are  shown.  Eddy 
currents  of  very  high  density  exist  under  the  insulator.  Of  particular 
interest  is  the  circulation  component 

<»A  -  - 


From  Fig.  12  it  can  be  seen  that  large  values  of  Jz  exist  close  to  the  plane 
z  =  0  which  must  fall  to  zero  at  z  =  0  leading  to  large  values  of  near 

3z 

z  =  0.  The  eddy  currents  with  positive  circulation  (Fig.  20)  may  be 

associated  with  the  negative  J  present  in  the  same  region  (Fig.  12). 

z 


3.1.6  Conclusions 


An  investigation  of  the  three-dimensional  electrical  conduction 
field  in  a  slanted  electrode  wall  MHD  generator  has  shown  that  a  condition 
for  minimum  Joule  dissipation  on  the  insulator  section  of  the  slanted  wall 
is  that  the  tangential  electric  field  E£  be  made  to  vanish.  For  a  finely 
segmented  channel  (t/H  *  1/10) it  was  shown  that  a  tangential  equipotential 
co»il.d  be  imposed  on  the  insulating  section  whereas  for  a  coarser  segmentation 
(t/H  =  1/2)  the  tangential  equipotential  was  not  maintained.  In  the  case 
of  thin  velocity  boundary  layers,  it  was  shown  that  a  parasitic  internal 
circulating  current  exists  due  to  the  shorting  of  the  field  E'  by  the  conductors 
inclined  to  the  magnetic  field.  This  loop  current  which  dissipates  energy 
within  the  channel  without  producing  any  useful  power  output  degrades  the 
performance  of  the  generator.  This  performance  degradation  is  evident  in 
the  comparison  of  the  internal  imped  mce  Hi/(Ri  )ideal  for  'thick  (1.19)  and 
tnin  (1.79)  boundary  layers.  The  effect  of  the  parasitic  current  may  be 
diminished  by  selecting  geometries  with  H  >>  a. 


The  original  assumption  that  the  slanted  electrode  wall  would 
relieve  the  high  current  densities  at  the  upstream  edge  of  the  electrode 
by  allowing  the  current  to  enter  the  electrode  along  the  magnetic  field 
lines  is  justified  by  the  detailed  calculation  of  the  current  distribution. 

In  Fig.  21  the  current  distribution  in  the  slanted  wall  geometry  at  the 
central  section  of  the  electrode  is  compared  with  results  for  parallel 
wall  channels.  Tt  can  be  seen  that  although  the  slanted  wall  distributions 
are  for  larger  Hall  parameters  3  =  It  they  are  more  uniform  than  parallel 
wall  distributions  at  lower  Hall  parameters  6  =  2,  3. 

A  distinguishing  feature  of  the  three-dimensional  generator  was 
the  existence  of  circulation  within  the  current  field  even  with  uniform 
conductivity.  This  circulation  was  shown  to  be  proportional  to  the  gradients 
of  the  current  in  the  magnetic  field  direction.  The  eddy  currents  observed 
in  the  apex  region  where  large  gradients  in  Jz  exist  were  an  indication  of 
the  circulation  present.  These  eddy  currents  are  confined  to  the  apex 
within  a  region  of  the  order  of  the  period  length  i,  and  hence,  for  geo¬ 
metries  with  £/a  »  1,  the  eddy  currents  may  be  confined  to  a  locally  small 
region  of  the  wall.  For  a  uniform  gas  in  such  geometries  they  should 
therefore  not  seriously  affect  the  overall  dissipation  on  the  electrode 
wall.  Indeed,  the  low  values  of  internal  impedance  obtained  even  with  the 
eddy  currents  present  (but  with  uniform  o,  B)  are  indicative  of  their  minimal 
effect.  In  a  nonequi librium  gas,  however,  the  high  dissipation  in  the  apex 
region  generated  by  the  eddy  currents  may  be  sufficient  to  locally  elevate 
the  conductivity  and  short  the  plasma  rendering  the  actual  segmentation 
ineffective.  Alternate  slanted  wall  arrangements  are  therefore  desirable 
in  which  sharp  changes  of  the  boundary  geometry  in  the  magnetic  field  di¬ 
rection  are  avoided  while  the  electrode  wadi  preserves  a  finite  angle 
relative  to  the  magnetic  field. 


IV.  REFERENCES 


1.  Research  on  New  Concepts  in  Energy  Conversion,  Annual  Technical 
Progress  Report  No.  1,  November  15,  1969. 

2.  Research  on  New  Concepts  in  Energy  Conversion,  Quarterly  Technical 
Progress  Report  No.  22,  March  15,  1970. 

3.  Kerrebrock,  J.  L.,  and  Draper,  J.  S.,  "Nonequilibrium  MHD  Generators 
with  Molecular  Gases,"  AIAA  Paper  No.  70-41,  January  19-21,  1970. 

4.  Hurwitz,  K.,  Sutton,  G.,and  Tamor,  S.,  "Influence  of  Tensor  Conductivity 
on  Current  Distribution  in  an  MHD  Generator,"  J.  of  Appl.  Phys.  32^,  2, 
February  1961. 

5.  Kerrebrock,  J.  L.,  AIAA  Journal,  Vol.  4,  November  1965,  p.  1938. 

6.  Bernard,  J. ,  Ricateau,  P. ,  and  Zettwoog,  P. ,  "Effect  of  Ionization 
Saturation  and  Duct  Shape  on  Losses  Near  Walls , "  Proceedings  of  Inter¬ 
national  Symposium  on  MHD  Electrical  Power  Generator,  Salzburg,  Austria, 
July  4-8,  1966. 

7.  Lowenstein,  A.  I.,  "The  Performance  of  a  Slanted  Electrode  MHD  Generator," 
M.S.  Thesis,  Dept,  of  Aero.  &  Astro.,  Mass.  Inst,  of  Tech.,  June  1970. 

8.  Crown,  J.  C.,  "Analysis  of  Magnetohydrodynamic  Generators  having  Segmented 
Electrodes  and  Anisotropic  Conductivity,"  United  Aircraft  Corporation 
Report  R-1352-2,  1961. 

9.  Oliver,  D.  A.,  SUIPR  Report  No.  163,  Institute  for  Plasma  Research, 
Stanford  University,  June  1967. 


15 


VW  '»<■!  *3 wjwwwt; ; 


Section  II 

figures  la,  lb:  A  view  of  the  coaxial-preionizer  discharge,  taken 
from  the  exit  of  the  generator.  The  lover  figure  (la)  is  a  cross- 
sectional  schematic  of  the  preionizer.  The  top  view  shows  the  pre¬ 
ionizer  in  the  stratified  mode,  (with  spinning  coil),  as  evidenced 
by  the  visibility  of  the  cathode  support  bar,  which  is  upstream  of 
*:he  discharge,  at  the  sides  of  the  passage.  In  the  second  figure, 
(without  spinning  coil),  the  discharge  fills  the  entire  channel. 
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Section  II 

Figures  3a ,  3b :  Shadowgraph  pictures  of  the  upstream  (3a)  and  downstream  (3b) 
portions  of  the  channel,  showing  weak  oblique  shocks  emanating  from  the  nozzle 
and  comparatively  strong  normal  shocks  in  the  last  one  fourth  of  the  channel. 
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Electrode 

frame  IOOA 


Electrical  discharge  system  for  pulsed  molecular  gas  experiment 
including  6 7 u F ,  9KV  supply,  initiating  and  crowbaring  circuits. 


Section  II 

Figure  7 :  The  discharge  tube,  magnet  structure,  frame  i  mirror 
mounts  for  the  pulsed  molecular  gas  experiment. 
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Section  III 

Figure  3;  X-Y  cross  section  of  channel  showing  m-n  mesh. 
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a)  Region  displayed  In  slantwall  figures. 
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b)  Representation  of  current  vector  entering  electrode. 

Section  III 

Fig.  5  Further  explanation  for  figures  displaying  current  vectors 
on  slantwall. 
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a)  phi  distribution 
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b)  3  distribution 
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Section 

:  |,  j  ,  D  distributions  on  slantwall. 
ky=1'2,  Ix=0,  fine  clectrode  segmentation,  thin  boundary  layer. 
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Section  III 

Fig.  1J  Current  vector  pattern  in  y-z  cross  section. 

k  -1/2,  Ix=0,  fine  eJcctrode  segmentation,  thin  boundary  lay 
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irrent  vectors  in  the  x-y  plane. 

.=1/2,  I  =0,  fine  electrode  segmentation,  thin  boundary  layer 
See  Fig.1*  for  further  explanation.) 


a)  phi  distribution 


'  ,  ‘  i  'Vn  \  N  \  " 

t  t  s  N  '  V  '  \ 

t  t  \  n  '  v  4  -y\ 

t  l  H  N  V  V  '  '  Ju\ 

t  i  '  v  v  v  -  ,\ 

t  1  \  N  '  V  '  -  ,  [//.  *\ 

i  v  n  *  N  ’  T/  t  ' 

t  v  \  \  '  '  4  f/jL,  /  / 

t  V  N  V  V  '  *  ^tfl  t  ,  ’ 

!  I  \  N  '  '  -  '  J/L,/  1 
,  VI  N  *  s  '  *  TT7  /  .  t 

,  \  v  \  '  •  "JIlj  i 


'.'.VA'v* 

-A'AV.-/. 

\  V  N  V  '  ‘  '  A 

\  V  S  V  '  '  '  ^ 

.  V  V  1  V  *  '  '  * 

V  N  >  V  V  '  ' 

v  v  v  : '  v  v  ' 

V  i  S  '  '  '  ^ 

\  \  i  V  '  '  '  ^ 

'  V' 

v\>'NSVV’ 
\\\  N  v  V  >  ^ 

\  \  '  V  >  v  v  " 

\  x  s  N  v  '  ^ 

\;  v  \  V  v 


tM| 

~ri 

i  i 

~r~t  t 

i  l  n 

T1  I 

VS ' ! 
rV  j 

rv 
v V  * 

fj—i ; 

a  '  J 
ft  X-  » 


Current  vectors  on  slantwall. 

k  =1/2,  1  =0,  coarse  electrode  segmentation,  thin  boundary  layer 
(^See  Fig*5  for  further  explanation.) 
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Fig.  15  Current  vector  pattern  in  y-z  cross  section 

k  *1/2,  I  “0,  coarse  electrode  segmentation,  thin  boundary 
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Fig.  19 


b)  Joule  dissipation  distribution 

3n,  D  distributions  on  slontvall.  Apex  analysis. 
ly=l/2,  Ix=0,  fine  electrode  segmentation,  tMr  be  in d pry 
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Fig.  21  Current  distribution  on  Electrode. 
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